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A Brownian particle's random motions can be rectified by a periodic potential energy landscape 
that alternates between two states, even if both states are spatially symmetric. If the two states 
differ only by a discrete translation, the direction of the ratchet-driven current can be reversed by 
changing their relative durations. We experimentally demonstrate flux reversal in a symmetric two- 
state ratchet by tracking the motions of colloidal spheres moving through large arrays of discrete 
potential energy wells created with dynamic holographic optical tweezers. The model's simplicity 
and high degree of symmetry suggest possible applications in molecular-scale motors. 



Until fairly recently, random thermal fluctuations were 
considered impediments to inducing motion in systems 
such as motors. Fluctuations can be harnessed, however, 
through mechanisms such as stochastic resonance p) and 
thermal ratchets Q , as efficient transducers of input en- 
ergy into mechanical motion. Unlike conventional ma- 
chines, which battle noise, molecular-scale devices that 
exploit these processes actually requite thermal fluctua- 
tions to operate. 

This article focuses on thermal ratchets in which the 
random motions of Brownian particles are rectified by a 
time-varying potential energy landscape. Even when the 
landscape has no overall slope and thus exerts no aver- 
age force, directed motion still can result from the accu- 
mulation of coordinated impulses. Most thermal ratchet 
models break spatiotemporal symmetry by periodically 
translating, tilting or otherwise modulating a spatially 
asymmetric landscape 2]. Inducing a flux is almost in- 
evitable in such systems unless they satisfy conditions of 
spatiotemporal symmetry or supersymmetry Even 
a spatially symmetric landscape can induce a flux with 
appropriate driving 0, 0, 0, 0] • Unlike deterministic mo- 
tors, however, the direction of motion in these systems 
can depend sensitively on implementation details. 

We recently demonstrated a spatially symmetric three- 
state thermal ratchet for micrometer-scale colloidal par- 
ticles implemented with arrays of holographic optical 
tweezers, each of which constitutes a discrete potential 
energy well Q. Repeatedly displacing the array first by 
one third of a lattice constant and then by two thirds 
breaks spatiotemporal symmetry in a manner that in- 
duces a flux. Somewhat surprisingly, the direction of 
motion depends sensitively on the duration of the states 
relative to the time required for a particle to diffuse the 
inter- trap separation [2| . The induced flux therefore can 
be canceled or even reversed by varying the rate of cy- 
cling, rather than the direction. This approach builds 
upon the pioneering demonstration of unidirectional flux 
induced by a spatially asymmetric time-averaged opti- 
cal ratchet 0, and of reversible transitions driven 
by stochastic resonance in a dual-trap rocking ratchet 
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FIG. 1: One complete cycle of a spatially-symmetric two-state 
ratchet potential comprised of discrete potential wells. 



Here, we demonstrate flux induction and flux reversal 
in a symmetric two-state thermal ratchet implemented 
with dynamic holographic optical trap arrays |l2t llflj . 
The transport mechanism for this two-state ratchet is 
more subtle than our previous three-state model in that 
the direction of motion is not easily intuited from the pro- 
tocol. Its capacity for flux reversal in the absence of ex- 
ternal loading, by contrast, can be inferred immediately 
by considerations of spatiotemporal symmetry. This also 
differs from the three-state ratchet [7j and the rocking 
double- tweezer [ToL 111) in which flux reversal results from 
a finely tuned balance of parameters. 

Figure ^ schematically depicts how the two-state 
ratchet operates. Each state consists of a pattern of dis- 
crete optical traps, modeled here as Gaussian wells of 
width a and depth Vq, uniformly separated by a distance 
L ^> a. The first array of traps is extinguished after time 
Ti and replaced immediately with a second array, which 
is displaced from the first by L/3. The second pattern 
is extinguished after time Ti and replaced again by the 
first, thereby completing one cycle. 

If the potential wells in the second state overlap those 
in the first, then trapped particles are handed back and 
forth between neighboring traps as the states cycle, and 
no motion results. This also is qualitatively different from 
the three-state ratchet, which deterministically transfers 
particles forward under comparable conditions, in a pro- 
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cess known as optical peristalsis 0, El ■ The only way 
the symmetric two-state ratchet can induce motion is if 
trapped particles are released when the states change and 
then diffuse freely. 
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FIG. 2: (a) Displacement function f(t). 
tilting-ratchet driving force, F(t) = —f]f{t). 
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The motion of a Brownian particle in this system can 
be described with the one-dimensional Langevin equation 



r,x(t) = -V'(x(t)-f(t)) + at), 



(1) 



where r\ is the fluid's dynamic viscosity, V(x) is the po- 
tential energy landscape, V'(x) = dV(x)/dx is its deriva- 
tive, and £ (t) is a delta-correlated stochastic force repre- 
senting thermal noise. The potential energy landscape in 
our system is spatially periodic with period L, 



V{x + L) = V(x). 



(2) 



The time- varying displacement of the potential energy in 
our two-state ratchet is described by a periodic function 
f(t) with period T = T 1 +T 2 , which is plotted in Fig.^a). 

The equations describing this traveling potential 
ratchet can be recast into the form of a tilting ratchet, 
which ordinarily would be implemented by applying an 
oscillatory external force to objects on an otherwise fixed 
landscape. The appropriate coordinate transformation, 
y{t)=x{t)-f{t) 12, yields 



w(*) = -v'(y(t)) + F(t) + £(t), 



(3) 



where F(t) = —T]f(t) is the effective driving force. Be- 
cause f(t) has a vanishing mean, the average velocity 
of the original problem is the same as that of the trans- 
formed tilting ratchet (x) = {y} , where the angle brackets 
imply both an ensemble average and an average over a 
period T. 

Reimann has demonstrated 0, Q that a steady-state 
flux, (y) ^ 0, develops in any tilting ratchet that breaks 
both spatiotemporal symmetry, 

V(y) = V(-y), and - F(t) = F(t + T/2), (4) 

and also spatiotemporal supersymmetry, 

-V(y) = V(y + L/2), and - F(t + At) = F(-t), (5) 
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FIG. 3: Steady-state drift velocity as a function of the relative 
dwell time, T 2 /Ti, for PV = 2.75, L = 5.2 ^m, a = 0.65 fj,m, 
and various values of T/t. Transport is optimized under these 
conditions by running the ratchet at T/t = 0.193. 



for any At. No flux results if either of Eqs. (0J or JSJ is 
satisfied. 

The optical trapping potential depicted in Fig. ^ is 
symmetric but not supersymmetric. Provided that F(t) 
violates the symmetry condition in Eq. J3J, the ratchet 
must induce directed motion. Although F(t) is super- 
symmetric, as can be seen in Fig- EJb) , it is symmetric 
only when T\ — T%. Consequently, we expect a parti- 
cle current for T\ ^ T 2 . The zero crossing at T\ = T 2 
furthermore portends flux reversal on either side of the 
equality. 

We calculate the steady-state velocity for this system 
by solving the master equation associated with Eq. (JTJ 
7, 15]. The probability for a driven Brownian particle to 
drift from position xq to within dx of position x during 
the interval t, is given by the propagator 

P(x, t\x ,0) dx = L(x - t,] dt ' 8{x - x ) dx, (6) 



where the Liouville operator is 
L(x, t 
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and where is the thermal energy scale |l5| . The 
steady-state particle distribution p(x) is an eigenstate of 
the master equation 



p{x) 



P{x,T\x ,Q) p(x )dx , 



and the associated steady-state flux is 
x — Xq 



T 



■ p(xq) P{x, T\xq, 0) dx dxo- 



(8) 



(9) 



The natural length scale in this problem is L, the inter- 
trap spacing in either state. The natural time scale, r = 
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L 2 /(2D), is the time required for particles of diffusion 
constant D to diffuse this distance. 

Figure [3] shows how v varies with Ti /T 2 for various 
values of T/t for experimentally accessible values of Vo, 
a, and L. As anticipated, the net drift vanishes for T\ = 
T 2 . Less obviously, the induced flux is directed from each 
well in the longer-duration state toward the nearest well 
in the short-lived state. The flux falls off as 1/T in the 
limit of large T because the particles spend increasingly 
much of their time localized in traps. It also diminishes 
for short T because the particles cannot keep up with the 
landscape's evolution. In between, the range of fluxes can 
be tuned with T. 

We implemented this model for a sample of 1.53 /im 
diameter colloidal silica spheres (Bangs Laboratories, lot 
number 5328) dispersed in water, using potential en- 
ergy landscapes created from arrays of holographic op- 
tical traps pi |l2t |l3j, ■ The sample was enclosed in a 
hermetically sealed glass chamber roughly 40 /im thick 
created by bonding the edges of a coverslip to a micro- 
scope slide, and was allowed to equilibrate to room tem- 
perature (21 ± 1°C) on the stage of a Zeiss S100 2TV 
Axiovert inverted optical microscope. A 100 x NA 1.4 oil 
immersion SPlan Apo objective lens was used to focus 
the optical tweezer array into the sample and to image 
the spheres, whose motions were captured with an NEC 
TI 324A low noise monochrome CCD camera. The mi- 
crograph in Fig.^Ja) shows the focused light from a 5 x 20 
array of optical traps formed by a phase hologram pro- 
jected with a Hamamatsu X7550 spatial light modulator 
|l7| . The tweezers are arranged in twenty-trap mani- 
folds 37 /im long separated by L = 5.2 ^m. Each trap is 
powered by an estimated 2.5 ± 0.4 mW of laser light at 
532 nm. The particles, which appear in the bright-field 
micrograph in Fig.^Jb), are twice as dense as water and 
sediment to the lower glass surface, where they diffuse 
freely in the plane with a measured diffusion coefficient 
of D = 0.33 ±0.03 /im 2 /sec. This establishes the charac- 
teristic time scale for the system of r = 39.4 sec, which 
is quite reasonable for digital video microscopy studies. 
Out-of-plane fluctuations were minimized by focusing the 
traps at the spheres' equilibrium height above the wall 

a 

We projected two-state cycles of optical trapping pat- 
terns in which the manifolds in Fig. 0{a) were alter- 
nately displaced in the spheres' equilibrium plane by L/3, 
with the duration of the first state fixed at T% = 3 sec 
and T 2 ranging from 0.8 sec to 14.7 sec. To measure 
the flux induced by this cycling potential energy land- 
scape for one value of T 2 , we first gathered roughly two 
dozen particles in the middle row of traps in state 1, as 
shown in Fig. 0{b), and then projected up to one hun- 
dred periods of two-state cycles. The particles' motions 
were recorded as uncompressed digital video streams for 
analysis |l9| . Their time-resolved trajectories then were 
averaged over the transverse direction into the prob- 
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FIG. 4: (a) Image of 5 x 20 array of holographic optical traps 
at L — 5.2/im. (b) Video micrograph of colloidal silica spheres 
1.53 /jm in diameter trapped in the middle row of the array 
at the start of an experimental run. (c) and (d) Time evolu- 
tion of the measured probability density for finding particles 
at T2 = 0.8 sec and T2 = 8.6 sec, respectively, with Ti fixed 
at 3 sec. (e) Time evolution of the particles' mean position 
calculated from the distribution functions in (c) and (d). The 
slopes of linear fits provide estimates for the induced drift ve- 
locity, which can be compared with displacements calculated 
with Eq. CHJ for /3V = 2.75, and a = 0.65 fim. (f) Mea- 
sured drift speed as a function of relative dwell time T2/T1, 
compared with predictions of Eq. . 



ability density, p(x,t)Ax, for finding particles within 
Ax = 0.13 /im of position x after time t. We also tracked 
particles outside the trapping pattern to monitor their 
diffusion coefficients and to ensure the absence of drifts 
in the supporting fluid. Starting from this well-controlled 
initial condition resolves any uncertainties arising from 
the evolution of nominally random initial conditions [7j. 

Figures OJc) and (d) show the spatially-resolved time 
evolution of p{x,t) for T 2 = 0.8 sec < T\ and T 2 = 
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8.6 sec > T±. In both cases, the particles spend most of 
their time localized in traps, visible here as bright stripes, 
occasionally using the shorter-lived traps as springboards 
to neighboring wells in the longer-lived state. The mean 
particle position (x(t)) = J x p(x,t) dx advances as the 
particles make these jumps, with the associated results 
plotted in Fig.HJe). 

The speed with which an initially localized state, 
p(x,0) fts S(x), advances differs from the steady-state 
speed plotted, in Fig. but still can be calculated as 
the first moment of the propagator, 



x(t)) = / yP(y,t\0,0)dy. 



(10) 



Numerical analysis reveals a nearly constant mean speed 
that agrees quite closely with the steady-state speed from 
Eq. ©■ 

Fitting traces such as those in Fig.^Je) to linear trends 
provides estimates for the ratchet-induced flux, which 
are plotted in Fig. 0Jf). The solid curve in Fig. 0Jf) 
shows excellent agreement with predictions of Eq. (|1C)|> 
for PV = 2.75 ± 0.5 and a = 0.65 ± 0.05 |xm. 
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FIG. 5: Toy model of diffusive molecular motor. 

Our implementation of the two-state ratchet involves 
updating the optical intensity pattern to translate the 
physical landscape. However, the same principles can be 
applied to systems in which the landscape remains fixed 
and the object undergoes cyclic transitions between two 
states. Figure CD depicts a model for an active two-state 
walker on a fixed physical landscape that is inspired by 
the biologically relevant trans por t of single myosin head 
groups along actin filaments [20j. The walker consists 
of a head group that interacts with localized potential 
energy wells periodically distributed on the landscape. 
It also is attached to a lever arm that uses an external 
energy source to translate the head group by a distance 
somewhat smaller than the inter- well separation. The 
other end of the lever arm is connected to the payload, 
whose viscous drag would provide the leverage necessary 
to translate the head group between the extended and re- 
tracted states. Switching between the walker's two states 
is equivalent to the two-state translation of the poten- 
tial energy landscape in our experiments, and thus would 
have the effect of translating the walker in the direction of 
the shorter-lived state. A similar model in which a two- 
state walker traverses a spatially asymmetric potential 
energy landscape yields deterministic motion at higher 



efficiency than the present model |21j. It does not, how- 
ever, allow for reversibility. The length of the lever arm 
and the diffusivity of the motor's body and payload deter- 
mine the ratio T/r and thus the motor's efficiency. The 
two-state ratchet's direction does not depend on T/t, 
however, even under heavy loading. This differs from the 
three-state ratchet 0, in which T/t also controls the di- 
rection of motion. This protocol could be used in the 
design of mesoscopic motors based on synthetic macro- 
molecules or microelectromechanical systems (MEMS). 
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